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We present results for the charged-particle multiplicity distribution at midrapidity in Au - Au 
collisions at yfs NN — 130 GeV measured with the PHENIX detector at RHIC. For the 5% most 
central collisions we find dN c h / dn\ v=0 — 622± l(stat) ±41(syst). The results, analyzed as a function 
of centrality, show a steady rise of the particle density per participating nucleon with centrality. 

PACS numbers: 25.75.Dw 



The Relativistic Heavy-Ion Collider (RHIC) at 
Brookhaven National Laboratory started regular oper- 
ation in June 2000, opening new frontiers in the study of 
hadronic matter under unprecedented conditions of tem- 
perature and energy density. The research is focused on 
the phase transition associated with quark deconfincment 
and chiral symmetry restoration expected to take place 
under those conditions. 

In this letter, we report results for the charged-particle 
multiplicity distribution at midrapidity in Au - Au col- 
lisions at y/s NN — 130 GeV, as measured with the 
PHENIX detector. These are the first RHIC results to 
span a broad impact parameter range. 

Particle density at midrapidity is an essential global 



variable for the characterization of high energy nuclear 
collisions, providing information about the initial condi- 
tions, such as energy density. The results presented here 
should help to constrain the wide range of theoretical 
predictions Q available at RHIC energies and to discrim- 
inate among various mechanisms of entropy and particle 
production. In particular, we analyze the particle density 
as a function of centrality, expressed by the number of 
participating nucleons. Such an analysis may shed light 
on the relative importance of soft versus hard processes 
of particle production and test the assumption of gluon 
saturation expected at RHIC energies ||||. Our results 
are compared to different models, to similar studies ob- 
tained in Pb-Pb collisions at the CERN SPS and to 
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a recent measurement performed at RHIC by PHOBOS 

The PHENIX detector is decribed in ref. §]. The 
present analysis relies primarily on three PHENIX sub- 
systems: two layers of Pad Chambers (PC), called PCI 
and PC3, used to determine the charged particle mul- 
tiplicity, the Zero Degree Calorimeters (ZDC) and the 
Beam-Beam Counters (BBC), used to derive the trigger 
and the off-line event selection. The PC provide three- 
dimensional coordinates along the charged-particle tra- 
jectories |9j. The two layers are mounted at radial dis- 
tances of 2.49 m and 4.98 m, respectively, from the in- 
teraction region. Each layer has 8 wire chambers with 
cathode pad readout (see Fig. 1) and covers 90° in az- 
imuth (4>) and ±0.35 units of pseudorapidity (77). The 
ZDC are small transverse-area hadron calorimeters that 
measure neutron energy within a 2 mrad (|?7| > 6) cone 
around the beam direction and are located at ±18. 25m 
from the center of the interaction region jl0|]. The BBC 
comprise two arrays of 64 photomultiplier tubes each 
equipped with quartz Cherenkov radiators. The BBC 
are located around the beam direction at ± 1.44m from 
the center of the interaction region covering the full 2ir 
azimuth and the range rj = ± (3.0 - 3.9) ]rj| . 

The primary interaction trigger is generated by a co- 
incidence between the two BBC with at least two pho- 
tomultipliers fired in each of them and a requirement on 
the collision vertex position, usually \z\ < 20 cm. Based 
on detailed simulations of the BBC, this trigger reflects 
[92 ± 2(syst)]% of the nuclear interaction cross section of 
7.2 barns p2| . Another trigger is generated by a coin- 
cidence between the two ZDC, each one with an energy 
signal larger than 10 GeV. This trigger reflects the nu- 
clear interaction plus the mutual Coulomb dissociation 
cross sections. Most BBC triggers (97.8%) also satisfy 
the ZDC trigger requirement. The small percentage of 
exclusive BBC triggers is due to inefficiencies of the ZDC 
trigger and background interactions. These two sources 
are difficult to distinguish and we estimate a background 
event contamination of [1 ± l(syst))% of the total event 
sample. 

This analysis is based on a sample of 137,784 events 
taken without magnetic field, satisfying the BBC trigger 
and with a reconstructed vertex position \z\ < 17 cm. 

The number of primary charged particles per event 
is determined on a statistical basis by correlating hits 
in PCI and PC3. The analysis procedure was devel- 
oped and corroborated by extensive simulations using 
the GEANT Q response of the PHENIX detector to 
events generated with HIJING [|l4j . The vertex is recon- 
structed using the following algorithm: all hits in PC3 
are combined with all hits in PCI and the resulting lines 
are projected onto a plane through the beam line, per- 
pendicular to the symmetry axis of the chambers (see 
Fig. 1). For events with more than ^5 tracks, the dis- 
tribution of these projections along the Z axis produces 



a distinct peak which defines the vertex position. For 
low-multiplicity events the vertex is reconstructed from 
the time difference between the two BBC. 

Once the vertex is known, all hits in PC3 are again 
combined with all hits in PCI and the resulting tracks 
are projected onto the plane previously defined. The dis- 
tribution of the distance R of the intersection points to 
the vertex position, is shown in Fig. 2. This distribution 
contains real tracks and tracks from the obvious combi- 
natorial background inherent to the adopted procedure 
of combining all hits in PC3 with all hits in PCI. The 
latter can be determined by a mixed event technique; in 
the present analysis, each sector in PCI was exchanged 
with its neighbor and the resulting combinatorial back- 
ground is shown in Fig. 2 by the dotted line. The yield of 
this background increases quadratically with R (leading 
to a linear dependence in the differential dN/dR vs. R 
presentation of Fig. 2). 

The R distribution of real tracks obtained by subtract- 
ing the background (Ns) from the total number of tracks 
is shown in Fig. 2 by the dashed line. The sharp peak 
at small R is due to tracks from primary particles orig- 
inating at the vertex and the long tail is due to decay 
products of primary particles decaying in flight. In prac- 
tice, the track counting is performed up to a given R 
value. The fraction of counted tracks as a function of R 
is obtained from the integral of the dashed curve normal- 
ized to the total integral up to R = 00. The tail at large 
R values is very well described by an exponential func- 
tion and therefore the extrapolation of the dashed curve 
to R = 00 is straightforward. The larger the value of R, 
the smaller is the correction for the fraction of uncounted 
tracks but the larger the background to be subtracted. 
Since the background can be reliably subtracted we per- 
formed the track counting up to R=25cm, thus including 
95.9% of all tracks. 

After subtracting the background, a raw multiplicity 
distribution is obtained to which we apply several cor- 
rections to obtain the corrected distribution of primary 
charged particles: 

(i) A correction of 15.3% accounts for inactive gaps 
between the chambers, inactive electronic readout cards 
and dead pads in the PCI and PC3 detectors. 

(ii) Using cosmic rays, the pad chamber hit efficiency 
was measured to be 99.4% for isolated single hits in agree- 
ment with an analytical study of the chamber perfor- 
mance |fL5[ . 

(iii) The track losses due to the finite double hit reso- 
lution of the chambers depend on the event multiplicity. 
The losses occur in the direct counting of tracks and in 
the combinatorial background subtraction. The two ef- 
fects were studied in great detail with Monte Carlo tech- 
niques. For the first one, the correction amounts to 13.3% 
for the 5% most central collisions, whereas it is only 4.9% 
for the 20-25% bin. The second effect is a reduction of 
the combinatorial background which leads to an increase 
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in the track multiplicity by 0.036 x Nb- 

(iv) Finally there is a correction due to uncounted 
charged tracks. Two sources contribute here. First, there 
is a correction of 4.3% from tracks missed because the 
analysis is limited to R=25cm. The second source has 
two components. On the one hand there are primary 
charged particles (mainly n^) which decay in flight. A 
large fraction of these decays is accounted for since, as 
discussed above, they produce the tail of the R distri- 
bution in Fig. 2. However, there is still a small number 
of decays which miss altogether PCI and/or PC3 and 
those have to be added. On the other hand, there is 
feed-down from neutral particle decays (mainly K° s and 
7r°) which lead to valid tracks. Those have to be sub- 
tracted. Both components depend on particle compo- 
sition and momentum distribution. Lacking precise in- 
formation about them, we performed Monte Carlo sim- 
ulations, applying to the GEANT response of the pad 
chambers the same analysis procedure as to real data 
and comparing the resulting multiplicity to the original 
HIJING input. The result is a net correction of only 
~2.8%. This correction is quite robust against changes 
in the input. For example, a uniform 20% increase in 
momentum p of all pions reduces the net correction to 
1%. 

The systematic error in the multiplicity due to correc- 
tions (i) and (ii) is estimated at less than 2%. To these 
we have to add the errors associated with the background 
subtraction, double hit resolution (which are multiplicity 
dependent) and particle decays (which is multiplicity in- 
dependent). We estimate the first two uncertainties to be 
4.6% at the highest multiplicities, based on Monte Carlo 
guidance and a comparison with an identical analysis in 
which we imposed a 50% larger double hit resolution in 
PCI and PC3. For the particle decay correction, we as- 
sign an error of 4% based on the effect of varying the 
momentum distribution and the particle composition in 
the simulations and a comparison with the results ob- 
tained with another data set measured while the detec- 
tors were retracted by 44 cm from the nominal position. 
Adding these errors in quadrature results in a total sys- 
tematic error of 6.5% at the highest multiplicities. An 
additional error in dNch/df], which enters in the analysis 
of the multiplicity versus centrality presented below, is 
the uncertainty in the total number of events due to the 
BBC trigger (±2%) and the event contamination (±1%). 

After applying the appropriate corrections we obtain 
in the lower panel of Fig. 3 the minimum-bias charged- 
particle multiplicity distribution, in the track acceptance 
|r?| < 0.35, 5<j) = 88.4°. A factor of 5.82 thus converts the 
observed number of tracks to dN c h / drj\ n= Q in one unit of 
pseudorapidity and full azimuth, yielding the lower hor- 
izontal scale in Fig. 3. 

Figure 3 also shows the multiplicity distributions for 
the four most central bins, 0-5% to 15-20%. The bins 
were defined by cuts in the space of BBC versus ZDC 



analog response (see Fig. 3 - upper panel) and refer to 
percentiles of the total interaction cross section of 7.2 
barns. In order to avoid the ambiguities inherent in the 
BBC vs. ZDC distribution, we selected events with in- 
creasing centrality based upon the monotonic response of 
another PHENIX detector. The centroids of these events 
projected onto the ZDC-BBC space determine the cen- 
trality contour indicated by the solid line. The cuts are 
made perpendicular to this contour. Simulations of the 
BBC and ZDC response were used to account for the ef- 
fect of physics and detector fluctuations in the definition 
of these event classes and to relate them via a Glauber 
model Jl^ l to the number of participating nucleons N p 
and of binary collisions N c . The average charged particle 
density scaled up to one unit of rapidity and the corre- 
sponding numbers N p and N c are tabulated in Table 1 
for various centrality bins. 

The PHOBOS experiment has recently reported an av- 
erage dNch/dr/^Q = 555±12(stat)±35(syst) for the 6% 
most central collisions Q]. For the same centrality bin, 
we&addNch/driin^o = 609±l(stat)±37(syst). From the 
preliminary results reported by NA49 in Pb-Pb collisions 
at v^njv — 17.2 GeV ||, we derive a particle density 
dN c h/ dy\ y = Q = 410. Our result, dN^/dr]^^ = 622 for 
the same centrality bin (0-5%), represents an increase of 
~ 80% after scaling it by 1.2 to account for the transfor- 
mation from rj to y. 

It has recently been emphasized [0 that the centrality 
dependence of dN c h/drj allows one to discriminate be- 
tween various models of particle production. We show in 
Fig. 4 our results for dN c h/dr) per participant pair as a 
function of the number of participants N p . Fig. 4 also 
shows the pp value at the same y/s taken from the UA5 
analysis |[(|. It is interesting to note that the extrapo- 
lation of our data points to low multiplicities approaches 
the pp value. 

Models such as HIJING predict that there is a com- 
ponent of particle production from soft interactions that 
scales linearly with N p and a second component from 
hard processes (pQCD jets) that scales with N c . Fol- 
lowing that, we fit the data of Fig. 4 with the function 
dN ch /drj = A x N p + B x N c using the values of N c and 
N p tabulated in Table 1. The results of the fit are shown 
in Fig. 4. Note that the errors in A and B are anticorre- 
lated. In such models, the values of A and B imply a large 
contribution of hard processes to particle production, 
which increases with centrality from ~ 30% at 7V P =68 to 
~ 50% for central collisions. HIJING predicts the same 
trend although the calculated values are lower than the 
data by ~ 15%. However, this is not a unique inter- 
pretation. At Alternating Gradient Synchrotron (AGS) 
energies, where hard processes do not occur, the parti- 
cle production per participant was also observed to in- 
crease with centrality 0. At the CERN SPS energy of 
^s NN = 17.2 GeV, a similar behavior was observed. Us- 
ing a different parametrization, dN c h/d-q cx N p , WA98 
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finds a best fit value of a = 1.07 ± 0.04 [|. Experi- 
ment WA97 quotes a value a = 1.05 ± 0.05, consistent 
with the result of WA98 but also compatible with their 
assumption of proportionality between multiplicity and 
participants ||. A good fit to our data can also be ob- 
tained with this functional form with a higher value of 
a = 1.16 ±0.04. One should note that the CERN results 
are in the lab frame whereas our are in the center of mass 
system. 

Other models such as the EKRT || predict that at 
RHIC energies, the large production of semi-hard gluons 
in a small volume may saturate the gluon density. The 
resulting gluon fusion limits the total entropy production 
and thus lowers the final particle production per partici- 
pant. The predictions of the EKRT model are also shown 
in Fig. 4. We observe no such saturation effect within our 
errors for Au + Au collisions at ^/s^ = 130 GeV, but 
instead see a steady rise in the particle production per 
participant pair. 
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FIG. 1. Sketch of the geometry. For clarity, three PC3 
sectors have been removed from the drawing. 
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FIG. 2. Tracks per event and per cm as a function of the 
distance to the event vertex. 
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FIG. 3. BBC vs. ZDC analog response (top panel) and 
minimum-bias multiplicity distribution in the PHENIX mea- 
surement aperture (lower panel). The lower axis converts the 
observed distribution to the corresponding average dN c h/dr) 
for track multiplicities less than ~ 120; beyond that value the 
shape of the distribution has a significant contribution from 
fluctuations into the measurement aperture. 
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FIG. 4. Charged-particle pseudorapidity density per par- 
ticipant pair vs. the number of participants. Predictions 
from HIJING § and EKRT § models, and a simple phe- 
nomenological fit are also shown. The shaded area represents 
the systematic errors of dN c h/dri and N p . The errors of N c 
are in Table 1. 



TABLE I. Charged-particle density, number of partici- 
pants and collisions for various centrality bins expressed as 
percentiles of a geo —7.2 b. The statistical errors are negligible 
and only systematic errors are quoted. 



Centrality bin 


< dN ch /dr)i v=0 > 


< N p > 


<N C > 


- 5 


622 ± 41 


347 ± 10 


946 ± 146 


5 - 10 


498 ± 31 


293 ± 9 


749 ±116 


10 - 15 


413 ± 25 


248 ± 8 


596 ± 93 


15 - 20 


344 ± 21 


211 ± 7 


478 ± 75 


20 - 25 


287 ± 18 


177 ± 7 


377 ± 61 


25 - 30 


235 ± 16 


146 ± 6 


290 ± 47 


30 - 35 


188 ± 14 


122 ± 5 


226 ± 38 


35 - 40 


147 ± 12 


99 ± 5 


170 ± 30 


40 - 45 


115 ± 11 


82 ± 5 


130 ± 24 


45 - 50 


89 ± 9 


68 ± 4 


101 ± 19 
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